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Abstract. 

In the theory of resonant scattering, the double differential cross section involves the 
computation of a multifold integral of a 4-point correlation function, which generalizes 
the traditional 2-point correlation function of Van-Hove for potential scattering. In 
the case of a neutron-crystal interaction, the numerical computation of these multifold 
integrals is cumbersome. In this paper, a new approximation is suggested. It is 
based on a factorization of the differential cross section into one function describing 
the exchange of kinetic energy between the neutron and the bound nucleus (phonons 
dynamic) and a function related to the nuclear scattering amplitude. This formalism 
is then applied to the modeling of resonant scattering of a neutron by ^ssjj jj-^ ,^ UO2 
crystal lattice. 
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1. Introduction 



The computation of the double differential scattering cross section (DDCS) at low 
incident neutron energies is required to solve neutron transport problems. Wigner and 
Wilkins [1] used a two-body kinematic approach to study potential scattering in a free 
gas. Under the same classic assumptions, Blackshaw and Murray [2] studied the case 
of an energy- dependent cross section and further generalizations were investigated by 
Ouisloumen and Sanchez and Rothenstein and Dagan [1]. 

A general quantum formalism due primarily to Van Hove [5] expresses the DDCS of 
potential scattering by a bound nucleus as a Fourier transform of a 2-point correlation 
function. Kazarnovski et al. ^ and Word and Trammell ^ extended the Van Hove 
theory to study resonant processes. The resonant DDCS becomes a Fourier double- 
Laplace transform of a 4-point correlation function. 

In the case of a harmonic crystal, the numerical computation of the multifold 
integral is difficult because of the highly oscillating behavior of the 4-point correlation 
function. A simplified multiphonon model, known as the "uncoupled phonon 
approximation" (UPA) was proposed by Naberejnev [9J. However, the validity of 
this approximation is questioned especially within the short time range (case of high 
temperature for instance). The present work suggests a new approximation that gives 
a correct limit at short time and provides a simple formula to compute the scattering 
kernel. 



2. General formalism 



2.1. The 4-point correlation function 



When the scattering amplitude is independent of incident neutron energy, it is known 
since the work of Van Hove [5j that the scattering DDCS is related to the Fourier 
transform of a two point correlation function X2{''~) '■ 



(Pa <^p kf 



dVtdEf Arc ki 



^to(r)}(AE). (1) 



AE = Ef — Ei is the energy transfer with Ei and Ef the energy of the neutron 
before and after the scattering event. The Fourier transform is defined as 

HX2{r)}{x) = e—X2{r)dT. (2) 

^{X2{t)}{AE) is the usual scattering function S{a,(3) = S{Ak, AE) which 
depends on the momentum transfer Ak = kj — ki where ki and kj are the initial and 
final neutron wave vector respectively . The theory was later generalized to treat the 
case of resonant scattering [7]. In the incoherent approximation, the scattering cross 
section includes the resonant, the potential- resonant interference and potential terms. 
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It can be demonstrated that the DDCS of the resonant term for instance is related to 
the Fourier-double-Laplace transform £2-^ of a 4-point correlation function X4(r, 

E'o, r„ and F are the energy, the neutron and total width of the resonance. 
am = '^T^g^n/kf^ and g is the usual statistical spin factor. The Fourier-double-Laplace 
transform is defined as 

1 /"OO /"CX) /"OO 

C2J^{f{T,t,t')}{x,zl,z2) = — / / e-'^'^e-'^'e-''' f{T,t,t')dTdtdt'. (4) 

2n J-00 JO Jo 

Similarly, the resonant-potential interference term is expressed as a Fourier-single- 
Laplace transform of a 3 point correlation function X3(r, t). 

A quantum calculation shows that X4.{^j't,t') is a function of time- dependant 
displacement operators: 



X4(r, t, t') = < exp (^—iki.f{T + t — t')^ exp (ikf.f{T + t)^ 

X exp (^—ikf.f{t)^ exp (j.ki.f{0)^ >t ■ (5) 

Note that X2(t) = Xii'Ty 0, 0) and Xsi'^^ t) = Xii^^ ^) 0). < X >t denotes the thermal 
average at temperature T. In this paper, time r, t and t' are expressed in unit of energy 
and h = 1. The right hand side of Equation 1 and 3 should be multiplied by ( ^^^"^ ) 
which is implicitly assumed in the paper. M and m are the mass of the target nucleus 
and incident neutron and A = M/m respectively . The scattering is assumed isotropic 
in the center-of-mass frame and multiple scattering effects are neglected (single collision 
approximation). Equation 3 assumes that the Hamiltonian of the target and compound 
nucleus is the same. This condition neglects the mass change of the target nucleus when 
the neutron is absorbed and re-emitted: M » m. 



2.2. Free-gas versus harmonic crystal-lattice model 

In the case where the target nucleus behaves like a free gas, the displacement operator 
is proportional to the momentum operator r{t) = r(0) +p{0)t/M. After averaging over 
a Maxwellian distribution of momentum at temperature T, xa within the free gas model 
(FGM) becomes 

X,{r,t,t') = exp{^[r(% - + {t - t')k^] - ^\-hit - t') + r(fc; - h^f}- (6) 

When the target nucleus is bound to a crystal lattice, the 4-point correlation 
function is more complex. Under the harmonic approximation, the Bloch theorem [H] 
transforms X4 into 
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X4(i-, t, t') = exp{-p{ki, ki, 0) - p{kf, kf, 0)} 

X exp{p(A^i, kf, -t') + p{kf, ki, +t) + p{kf, kf, r)} 

X exp{p(A;i, ki, t - t' + t) - p{ki, kf, t - t') - p{kf, k,, r + t)]. (7) 
For a cubic lattice with a phonon density of states pioj), 

k', +t) = — duj^-^\coth{uj /2kT)cos{ujt) - isin{ujt)]. (8) 
2M Jo OJ 



In the short time approximation for r, t, and t' the DDCS can be simphfied using 

2M 2M' 



p{k, k', +t) ^ p{k, k', 0) + k.k'i^^t' - -^), (9) 



1 f°° 

kTeff = 77 / dujp{u)ijiJCoth{uj /2kT). (10) 



with 

2 Jo 

Plugging [9] into [TJ we get the free gas formula [6] with T^ff instead of T. 
Consequently, the crystal-lattice model lead to the free gas model, with an effective 
temperature, when r, t, and t' are small (short time approximation). This is the same 
effective temperature as Lamb [I2] derived for the Doppler broadening of a capture 
resonance in a solid in the weak binding limit. 

2.3. The uncoupled phonon approximation 

The UPA approximation proposed in [9j was an attempt to compute the resonant 
scattering kernel for a harmonic crystal. It neglected t and t' in the coupling terms 
p{ki,ki,T — t' + t), p{ki,kf,T — t') and p{kf,ki,T + t) in Equation [7] and applied a 
short time approximation to p(ki, kf, —t') and p{kf, ki, +t) 0. The correlation function 
becomes 



xV\t, t,t')= exp (p{£k, t) - p(Ak, 0) 



The integral over r and t, t' can be performed separately. Since the first exponential 
is the Van- Hove function X2(t), the differential cross section can be factored out as 



d^a 1 k 



f 



dVtdE 4:71 ki 



X S{Ak, AE) X a{ki, kf). (12) 



The detailed mathematical expression for the so-called UPA cross section a{ki, kf) 
can be found in [9]. The UPA model separates the phonon dynamic and the nuclear 

I A refinement of the model called MUPA (modified uncoupled phonon approximation) computed 
p{ki, kf, —t') and p{kf, ki, +t) with a discrete phonon spectrum 
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interaction and provides a rather simple way to compute DDCS. However, within the 
UPA approximation, the 4-point correlation for small r, t, and t' becomes 

Xr^(r,t,t') = exp{^[r(fc^-fc.)^ + (t-t')^.^/]-^[fc.^/(t^+t'^)+r^(^/-fc.)^]}, (13) 

which differs markedly from the required equation [6] discussed in the previous 
section. Consequently, the DDCS within the UPA approximation does not give the 
correct limit for small t and t'. 



3. Proposed Model 

Our model seeks to get a factorization of the DDCS similar to the UPA model: 

"^'"^ -^^ X SiAk, AE) X aih, kj). (14) 



dQdE 47r k 

However, in order to get the correct limit at short time, we computed the a term 
using the short time approximation of the correlation function (Equation [6]), without 
the previous UPA approximation. 

For values of Ak ^ 0, the following notation is used. 



A 



'J'm^ .^2lE.-E.^iE,^E,'m-'§lUr. (15, 
\ ^ 2Ak Ak 2M 



It is found that the differential cross section, including resonant, potential and 
resonant-potential interference can be calculated almost analytically and can be split 
into nuclear and transfer terms. The details of the demonstration are presented in the 
Appendix. The final result is: 



d'^a _ I kf Qt/Ci r i r„ , , / r„ 



X S{Ak,AE) X [ (Xp + cr^—^p + d (Tpa„^—x ]■ (16) 



dndE Airh ^ ' ' , , ^ . , ^ 

For A ^ and Ak ^ 0, 



0r r xP r ^ P xP P 

ip = - — X— JCi— -,— -) y = X — L(— -, — -). (17) 
2 A ^2A 2A^ A ^2A 2A^ ^ ^ 



For A = and Ak ^ 



1 2x 

^ = ^ X 



1 + 1 + x^ 

where K{x,y) and L{x,y), are the real and imaginary part of complex 
complementary error function w{z) = exp[— (x + iyY] erfc [—i{x + iy)], related to 
the classic Voigt functions. 

K(.,y)^UH.^]=y-j-^0^at (19) 
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When Ak = 0, the same calculations lead to Equation [TS] with : 



A = x = 2\E,-E^-^]/V. (21) 

We recognize the usual Doppler-broadened cross section o"^, so that the differential 
cross section becomes 

- ^ X S{Ak,AE) X a^\E*), (22) 



where T* and E* depend on the cosine of the scattering angle /i = cos 9: 
kj{l-fi^) ^.hAk kikf 



T*{f^)=T f E*=E,-{Ef-E,)^-^, (23) 
and for Ak = 



Ak' Ak' 2M 



T* =T E* = Ei (24) 

Equation [23] is valid only when Ei — {Ef — Ei) — is positive. The previous 



equation can be further simplified noting that T* ^ ^'[(1 + /i)/2][l + AE / {2Ei)] ^ 
T(l + /i)/2 and when 1 — yU >> |AA;|/|A;j| we have ~ ^1/2 and this approximation 
produces satisfactory results when used over the whole range of values of /i. Making 
these approximations, we get a simple formula 



(Pa 1 k 



Equation [25] provides a simple way to compute the differential cross section within 
the free gas model. It also gives an approximate way to account for solid state effects 
by using the known scattering function S{Ak, AE) for the harmonic crystal. Note that 
Equation [25] was demonstrated for a single resonance and it is not known if this formula 
can be generalized to any form of free scattering amplitude. 



4. Application to neutron resonant scattering in UO2 

Low-enriched uranium oxide UO2 is widely used as nuclear-reactor fuel. With the 
present model, the scattering kernel of ^ss-q UO2 has been calculated near the first 
resonance at 6.67 eV by numerical integration of the differential cross section over 
scattering angles. The weighted phonon spectrum for ^ss-q UO2 has been taken 
from the measurement of Dolling et al. ^lOj. The spectrum features two acoustic modes 
around 14 and 21 meV. 

The S{Ak, AE) Van-Hove scattering function has been computed using the usual 
phonon expansion methods. The ^ss-q resonance parameters evaluated by Moxon and 
Sowerby [H] were used. 

The present model has been compared with the classic free gas kernel published by 
Blackshaw and Murray p] and studied by Ouisloumen and Sanchez [3]. Fig. 1 and 2 
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show the resulting scattering kernel at 300 K at two incident neutron energies, 6.674 
eV (the peak of the first resonance) and 6.520 eV which is the minimum of the cross 
section due the interference between the resonant and potential amplitudes. Compared 
to the UPA model (see Figures in Reference [9]), the present crystal model gives results 
closer to the free gas model. The only distinctive features are the small peaks around the 
elastic peak, associated with the one-phonon excitation of the acoustic modes. Note that 
when the S{Ak, AE) of the FGM is used in Equation [161 the numerical computations 
give the same results as Ouisloumen and Sanchez. 



5000 r 

^ 4000 




Secondary neutron energy (eV) 

Figure 1. Comparison of scattering kernels calculated with the present formalism 
(black line) and the free gas model of Ouisloumen and Sanchez [3] (grey line) at incident 
neutron energy of 6.674 eV. The dotted line is the kernel with a static model. 
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Secondary neutron energy (eV) 



Figure 2. Comparison of scattering kernels calculated with the present formalism 
(black line) and the free gas model of Ouisloumen and Sanchez [5] (grey line) at incident 
neutron energy of 6.520 eV. The dotted line is the kernel with a static model and the 
dashed line is the scattering cross section (cts/IOO in barns). 
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5. Conclusion 

An approximate formula is proposed to treat solid state effects in neutron-crystal 
interactions. The DDCS is the product of a Doppler-broadened scattering cross section 
and the usual Van-Hove scattering function. 

- — ^ X SiAk, AE) X {[e, + £^/]/2 - E,fi/A) . (26) 



d^ldE 4:71 ki 

This formula gives the correct free-gas limit for short-time range (high-temperature 
cases). A rigorous calculation of the DDCS within the crystal lattice model is still an 
open issue. The present work develops the case of a single isolated resonance but further 
generalizations to multilevel forms of collision matrix may be possible. 

In reactor applications, the model predicts small solid-state effects in the scattering 
of neutrons in UO2 at 300K, contrary to previous studies using the UPA approximation. 
Measurements of the secondary spectrum of neutrons scattered elastically in resonances 
would be valuable to check existing scattering models. 
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Appendix A. Derivation of Equation | 

For Ak 7^ 0, X4 in Equation [6] can be transformed into 



X4 = exp{^[( r - h^^_Jl I + (t - t%kj] 



2 ... 



kTj^ hAk{t-t')\ kfkj-ikjifj 



2M' 



2 , Ak'- ' ^ it-m}. (A.i) 

AA; / Ak 



^ ^ ^2 

Using the variable u = r — kiAk(t — t')/Ak , we recognize X2 the pair correlation 
function for a free gas: 

X2iu) = exr){—uAk^ - —u^Alz}. (A.2) 

Therefore, the Fourier transform of X4 (integration over r) will lead to the 
S{Ak, AE) terms with a phase factor to account for the change of variable t ^ u. 
The DDCS takes the form of Equation [22] with 

"P "P2 /.oo foo 

a = a^^— / e-''e-'*' (f){t-t')dtdt', (A.3) 

r 4 JO Jo 
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with z = -z [e, -Eo- {Ef - E,) ^ - + f and 



To compute the double Laplace transform, the following identity is used: 




(A.4) 



poo roo 

Jo Jo 



(f){t-t')dtdt' 



1 



roc roo 

/ e-''(p{t)dt+ / e-'''(p*{t)dt.{A.5) 



z + z* Jo 




Equation IA.3I becomes 



( 



2M 



kT 




dt 



4 Jo 



+ c.c. 



(A.6) 



c.c denotes the conjugate complex. This equation is then reduced into a single 
integral which can be further simplified using 



We recognize in I A . 6 1 and I A . 71 the real of the complex complementary error function. 
Equations [iHl [T7] and [20] are obtained in a similar way when the potential and potential- 
resonant interference terms are accounted for. 
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